Frequency-dependent counting statistics in interacting nanoscale conductors 
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We present a formalism to calculate finite-frequency current correlations in interacting nanoscopic 
conductors. We work within the n-resolved density matrix approach and obtain a multi-time cu- 
mulant generating function that provides the fluctuation statistics solely from the spectral decom- 
position of the Liouvillian. We apply the method to the frequency-dependent third cumulant of 
the current through a single resonant level and through a double quantum dot. Our results, which 
show that deviations from Poissonian behaviour strongly depend on frequency, demonstrate the 
importance of finite-frequency higher-order cumulants in fully characterizing transport. 

PACS numbers: 73.23.Hk,72.70.+m,02.50.-r,03.65.Yz 



Following the considerable success of shot-noise in the 
understanding of transport through mcsoscopic systems 
attention is now turning towards the higher-order 
statistics of electron current. The so-called Full Counting 
Statistics (FCS) of electron transport yields all moments 
(or cumulants) of the probability distribution P(n, t) of 
the number of transferred electrons during time t. De- 
spite their difficulty, measurements of the third moment 
of voltage fluctuations have been made [2 Jjl, and recent 
developments in single electron detection JUIjQ promise 
to open new horizons on the experimental side. 

The theory of FCS is now well established in the zero- 
frequency limit 0, H However, this is by no means 
the full picture, since the higher-order current correlators 
at finite frequencies contain much more information than 
their zero- frequency counterparts. Already at second or- 
der (shot-noise), one can extract valuable information 
about transport time scales and correlations. When the 
conductor has various intrinsic time scales like, for exam- 
ple, the charge relaxation time and the dwelling time of a 
chaotic cavity [lfj| . one needs to go beyond second-order 
in order to fully characterize electronic transport. Apart 
from this example, and some other notable exceptions 
Til . 13j |. the behaviour of finite-frequency correlators 
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beyond shot-noise is still largely unexplored. 

In this Rapid Communication, we develop a theory 
of frequency-dependent current correlators of arbitrary 
order in the context of the n-rcsolved density matrix 
(DM) approach, — a Quantum Optics technique 
that has recently found application in mcsoscopic trans- 
port [l5| . Within this approach, the DM of the sys- 
tem, p(t) , is unravelled into components p( n > (t) in which 
?i = n(t) = 0, 1, . . . electrons have been transferred to the 
collector. Considering a generic mesoscopic system with 
Hamiltonian H = TLs + Hl + ~Ht, where Hs and Hl 
refer to the system and leads respectively, and provided 
that the Born-Markov approximation with respect to the 
tunnelling term Ht is fulfilled, the time-evolution of this 



n-resolved DM can be written quite generally as 

p^{t) = C QP ^\t)+Cj P ^- 1 \t), (1) 

where the vector p^ n \t) contains the nonzero elements 
of the DM, written in a suitable many-body basis. The 
Liouvillian Cq describes the 'continuous' evolution of the 
system, whereas Cj describes the quantum jumps of the 
transfer of an electron to the collector. We make the infi- 
nite bias voltage approximation such that the transfer is 
unidirectional. By construction, this method is very pow- 
erful for studying interacting mesoscopic systems that are 
weakly coupled to the reservoirs, such as coupled quan- 
tum dots (QDs) in the Coulomb Blockade (CB) regime 
HI, IH or Cooper-pair boxes 17 1 . Within this frame- 
work, our theory of frequency-dependent FCS is of com- 
plete generality and therefore of wide applicability. 

In this n-resolved picture, electrons are transferred to 
the leads via quantum jumps and there exists no quan- 
tum coherence between states in the system and those in 
the leads. Thus, although the system itself may be quan- 
tum, the measured current may be considered a classical 
stochastic variable and therefore amenable to classical 
counting [3]. This observation allows us to derive vari- 
ous generalizations of classical stochastic results, and ob- 
tain a multi-time cumulant generating function in terms 
of local propagators. We illustrate our method by calcu- 
lating the frequency-resolved third cumulant (skewness) 
for two paradigms of CB mesoscopic transport: the sin- 
gle resonant level (SRL) and the double quantum dot 
(DQD). 

Equation ([1]) can be solved by Fourier transformation. 
Defining p(x,t) — J2 n (t)e mx , we obtain p(x,t) = 
M(x)p(x,t), with M(x) = C Q + e^Cj. Let N v be the 
dimension of M(x), and Aj(x); i = 1,...,N V , its eigen- 
values. In the x "~ > limit, one of these eigenvalues, Ai (x) 
say, tends to zero and the corresponding eigenvector gives 
the stationary DM for the system. This single eigen- 
value is sufficient to determine the zero-frequency FCS 
In contrast, here we need all N v eigenvalues. Using 
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the spectral decomposition, M(x) = V(x)A.(x)V -1 (x)j 
with A(x) the diagonal matrix of eigenvalues and V(x) 
the corresponding matrix of eigenvectors, the DM of the 
system at an arbitrary time t is given by 

p(x,t) = n(x,t-t )p(x,t ), (2) 

where fi(x; t) = e M M* = V(x)e Hx)t) V~ 1 (x) is the prop- 
agator in x-space, and p(x,to) is the (normalized) state 
of the system at to, at which time we assume no elec- 
trons have passed so that p^ n \to) = <5 n ,op(io) and thus 
p(x,to) = p(to)- The propagator in n-space, G(n,t) = 



-mx< 



£l(x,t), such that p {n) {t) = G(n,t - t )p(t ), 



for t > to, fulfils the property: G(n — no,t — to) = 
J2 n , G(n -n',t- t')G(ri - n , t' - t ), for t > t' > t , 
(n' = n(t'), n = n(i ))- This is an operator version of 
the Chapman-Kolmogorov equation [l9( . 

The joint probability of obtaining n\ electrons after 
t\ and n 2 electrons after t2, namely P > (ni, ti] n 2 , t 2 ) 
(the superscript '>' implies t 2 > ti), can be written in 
terms of these propagators by evolving the local probabil- 
ities P(n,t) = Tr{p(™)(£)} [201, and taking into account 
P(n 2 ,t 2 ) = J2 ni P > (ni,t 1 ;n 2 ,t 2 ), such that 

P > (n 1 ,t 1 ;n 2 ,t 2 ) = Tr {G{n 2 - m, t 2 - t%) 

xG(n 1 ,* 1 -t )p(to)}.(3) 

The total joint probability reads P(n\,t\,n 2 ,t 2 ) = 
TP>(n 1 ,t 1 ;n 2 ,t 2 ) = P > (n l ,t l ;n 2 ,t 2 )6(t 2 - h) + 
P < (n\,ti;n 2 ,t 2 )0(ti — t 2 ) where T is the time-ordering 
operator, 8(t) the unitstep function defined as 9(t) = 1 
for t > 0, and zero otherwise, and where P < (ni, t\\ n 2 , t 2 ) 
is the joint probability with t 2 < t\. It should be noted 
that, in contrast to the local probability P(n, t), the joint 
probability P(ni,ti;n 2 ,t 2 ) contains information about 
the correlations at different times. 

Result © may be alternatively derived using the Bayes 
formula for the conditional density operator [21| : 

P > (n 1 ,t 1 ;n 2 ,t 2 ) = P{n 1 ,t 1 )P > (n 2 ,f a |ni,ti) 

= Tr^ 



Tr 



{p^)(t 1 )}Tr{p^l" 1 )(i 2 )} 
{p(" 1 )(t 1 )}Tr|G(n 2 -n 1) t 2 -t 1 ) 



Tr 



{p(»0( tl )} 



The normalization in the denominator accounts for the 
collapse n = ni at t = t\ using Von Neumann's pro- 
jection postulate [2l|. Equation (J3J) is recovered when 



p( ni )(ti) is written as a time evolution from to- 

The two-time cumulant generating function (CGF) as- 
sociated with these joint probabilities is 

„-j r (xi,X2;ti,*2) = 



E 



P(n 1 ,t 1 ;n 2 ,t 2 )e iniXl+in2X2 , 



which, using Eq. ([3]), and e 



TTr(n( X 2,t 2 -t 1 ) 

xQ(xi + X2,ti-t )p{t )). (4) 



The above procedure can be easily generalized to obtain 
the N-timc CGF, which reads: 



N 



T Tr 



./c=l 



f/c! TjV-fc 



lp(*o) 



(5) 



where a k = Y,i=N+i-kXi, (xi,---,Xn), t 

:/'..... t n ) and = — tj, [22|. The multi-time CGF 
in Eq. ([5]) contains a product of local-time propagators, 
and expresses the Markovian character of the problem. 
It allows one to obtain all the frequency-dependent cu- 
mulants from the spectral decomposition of M(x). The 
iV-time current-cumulant (e = 1) is calculated using: 

S^(h,...,t N ) = (61(h)... SI(t N )) = 
= d tl ... dt N ((n(h) . ..n(t N ))) = 



= -(-i) N d tl 



■9t N d Xl 



d XN F(x;t] 



(6) 



The Fourier transform of with respect to the time 
intervals Tk, gives the Nth-order correlation functions 
as functions of N — 1 frequencies. In particular, the 
frequency-dependent skewness is a function of two fre- 
quencies which, as a consequence of time-symmctrization 
and the Markovian approximation, has the symmetries 
S^{uj,uj') = S^(uj',lj) = S( 3 > (w.w - w') = S^(lo' - 
u>,w') = S^(—w, —u)'), and is therefore real. The iVth- 
order Fano-factor is defined as FW = S^/(I). 

In the case where the jump matrix C j contains a single 
element, (Cj)ij = TftSiaSjp, which is the case for a wide 
class of models including our two examples below, all 
the correlation functions can be expressed solely in terms 
of the eigenvalues of M(0), and the N v coefficients 
c k = {V~ 1 C J V)kk = ^RV/3kV k ~J-. The second-order Fano 
factor then has the simple, general form 



F^(lj) = 1-2 



fc=2 



+ A 



(7) 



which has also been derived in other ways [23| . The skew- 
ness has the form F^{uj,uj') = -2 + E»=i f(2) (^) + 
F^(u),uj'), with i>i = oj, v 2 = to' , and ^3 = u — lo' . 
FW(u),w' ) is an irreducible contribution, the form of 
which is too cumbersome to be given here. The high- 
frequency limit of the skewness is F( 3 '(lj, 00) LO). 
As a first example we consider a SRL, described by p = 

t / — T L e ix Tfi \ 

(poO)Pn) and M(x) = ( -n n 1 1 m the basis 



Tl -Tr 

of 'empty' and 'populated' states, {|0), |1)}. Employing 
Eq. ([5|) , we obtain the known results for the current and 
noise, and arrive at our result for the skewness: 



Te * , gives F^(cj,cj') = 1-2T L T R 



uu ( r2 + -I) 

with 71 = Y\ + T%, 72 = 3r 2 , and T = T L + T R . The 
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FIG. 1: The third-order frequency-dependent Fano factor 
for the single resonant level, (a) Contour plot 
of F (3 \lo ,lu') as function of its arguments for Fa = Fl- 
(b) Sections F (3) (w,0) and F^ 1 (lo, -lo) show that the skew- 
ness is suppressed throughout frequency space both with 
respect to the Poissonian value of unity and to the noise 
F^ 2 '(u>). In contrast to the shotnoise, the skewness has a 
minimum at a finite frequency cj m , which exists in the cou- 
pling range (3 - VB) /2 < F L /F R < (3 + VB) /2. In di- 
rection lo = —lo', this minimum occurs, for Fl = Fr, at 
w m i/r H = VVio - 2/^3 whereas for u/ = the position 
of the minimum shifts slightly to LOmi/^R — 2/v3. (c) The 
maximum suppression occurs at Fl = Fr. 



zero-frequency limit F^ 3 ^(0, 0) agrees with Ref. j24j. 

This skewness is plotted in Fig. [T] from which the six- 
fold symmetry of is readily apparent. The third- 
order Fano factor gives, in accordance with the noise, a 
sub-Poissonian behaviour for all frequencies. This can 
be easily understood as a CB suppression of the long tail 
of the probability distribution of instantaneous current: 
due to the infinite bias, the distribution is bounded on the 
left by zero, but, in principle, is not bounded on the right 
(large, positive skewness). CB suppresses large current 
fluctuations which explains a sub-Poissonian skewness. 

Along the symmetry lines in frequency space with lo' = 
0, u) — lo and to = 0, the skewness is highly suppressed. 
In contrast with the noise, the minimum in the skewness 
occurs at finite frequency (Fig. QJ>) with the strongest 
suppression occuring at T L = Tr (Fig. [it). 

As a second example we consider a DQD in the strong 



CB regime [H, SfJ. In the basis of 'left' and 'right' 
states \L) and \R), which denote states with one excess 
electron with respect to the many body 'empty' state 
|0), the Hamiltonian reads H s = e(\L){L\ -\R)(R\) + 
T c (\L)(R\ + \R)(L\), with detuning e and coupling 
strength T c . The two levels \L), \R), arc coupled to their 
respective leads with rates Tl and T^. The DM vector 
is now p = (pooi Pll, PRR,^(p LR ),Im(p LR )) T , and the 
Liouvillian in this basis reads: 
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Comparison of the quantum-mechanical level-splitting 
A = + e 2 with the incoherent rates divides 

the dynamical behaviour of the system into two regimes. 
For A <C T^r, all eigenvalues of M(0) are real, and cor- 
respondingly, the noise and skewness are slowly-varying 
functions of their frequency arguments. In this regime, 
dephasing induced by the leads suppresses interdot co- 
herence, and the transport is largely incoherent. In the 
opposite regime, A 3> r^.^, two of the eigenvalues form 
a complex pair, A± w ±iA — r^/2 + 0(r/A) 2 , which sig- 
nals the persistence of coherent oscillations in the dots. 
The finite-frequency correlators then show resonant fea- 
tures at A, since these eigenvalues enter into the de- 
nominators, as in Eq. ([7]), giving rise to poles such as 
w=FA— iTr/2. The structure of the skewness is similar to 
the SRL for weak coupling (A <C Ti^j, but much richer 
in the strong coupling regime (A 3> ^l,r) (Fig.O- Now 
the skewness exhibits a series of rapid increases. From 
the origin outwards in the lu-uj' plane, we observe first a 
minimum at finite frequency and then inflexion points at 
lu ~ T R , \w\ = A, = A and \w - w'\ = A. Fig. 
shows sections in the lo-uj 1 plane, and the resonant be- 
haviour, in the form of Fano shapes, is most pronounced 
in F^> (u, —oS). Starting from high frequencies, the onset 
of antibunching occurs at u) = A. At higher frequencies, 
the system has no information about correlations and is 
Poissonian. The overall behaviour is seen in Fig. [^c) 
where we plot F^{uj 7 -lu) as a function of T c and w. 
The line lo = A delimits two regions: at high frequen- 
cies the skewness is Poissonian. At resonance, and after 
a small supcr-Poissonian region at lu > A, the system 
becomes sub-Poissonian (and even negative, for certain 
internal couplings). In the limit lo — > 0, our results qual- 
itatively agree with those of Ref. [13] for a noninteract- 
ing DQD: as a function of T c , the skewness presents two 
minima and a maximum (where the noise is minimum, 
not shown). In our case, however, the maximum occurs 
around A = Tr/2 — half that of the noninteracting case. 
Finally, we plot dF^(u>, —lo)/cIlo as a function of both e 
and lo (Fig.[5]i) where the resonances at w = A, lo = A/2 
and lo ~ Tr are clearly resolved. In contrast, the dcriva- 
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FIG. 2: Frequency-dependent Fano skewness for the dou- 
ble quantum dot in Coulomb blockade, (a) Contour plot in 
the strong coupling regime, T c = 3Fh, with r_L = Fr and 
e = 0. (b) Sections F (3) (w,0) and F {s) (uj, -w), and shot- 
noise .F^(w) show a series of abrupt increases with increas- 
ing uj. Both noise and skewness exhibit both sub- and super- 
Poissonian behaviour (c) Varying the internal coupling T c , 
the skewness shows rapid increases along the lines uj — A 
and uj = A/2. For u) > A the system is Poissonian (slightly 
super-Poissonian for uj > A), while for uj < A the transport is 
always sub-Poissonian. The skewness is strongly suppressed 
at low frequencies, (d) The derivative dF^ 3 '(uj,—uj)/doj as a 
function of frequency and detuning e for T c — 3Fl = 3r_R. 
Resonances occur at uj = A, A/2 and ~ Tr. 



tive dF^ (uj,0)/duj (not shown), exhibits a minimum at 
w = A for small e, which transforms into a maximum 
for e ^> T c . As expected, transport tends to be more 
Poissonian as e increases, signaling lost of coherence. 

Despite the simplicity of the models we have stud- 
ied, our results demonstrate the importance of finite- 
frequency studies. Deviations from Poissonian behaviour 
of higher-order cumulants are frequency-dependent, such 
that a comprehensive analysis in the frequency domain 
is needed in order to fully characterize correlations and 
statistics in electronic transport. 
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